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We consider an effective field theory for the nonleptonic decay in which a heavy quark 
decays into a pair of a heavy quark and antiquark having a small relative velocity and one 
relativistic (massless) quark. This effective theory is a combination of HQET, SCET, and 
a covariant modification of NRQCD. In the leading logarithm approximation the effective 
theory decay amplitude factorizes into the product of matrix elements of heavy-to-heavy 
and heavy-to-light currents. We discuss a possibility of factorization beyond the leading 
logarithm approximation and find it doubtful. The Wilson coefficients of the effective theory 
electro-weak (EWET) Lagrangian in the next-to-the leading logarithm approximation are 
calculated at the matching scale of the decay. The differential decay rate for the inclusive 
decay B — ► J/rfj + h in the effective theory framework is evaluated. 
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I. INTRODUCTION 

Decays of B mesons into charmonium are the subject of much theoretical and experimental 
study. For example, the CP asymmetry of the decay B — ► J/ipK s provides a very high precision 
determination of the unitarity angle (3. This is a clean determination because in the Standard 
Model (SM) of electroweak interactions the CP asymmetry does not depend on the decay rate. 
However, it is an interesting challenge in theoretical physics to determine from first principles the 
decay rate in the SM. Doubtless this would shed light into the underlying QCD dynamics. 

We propose to construct an effective field theory that captures the essential dynamics of decays 
of a B meson into charmonium plus light mesons. The theory we propose is a systematic expansion 
in small parameters and therefore the resulting approximations can be systematically improved. 
The building blocks of the effective theory are largely known: heavy quark effective theory (HQET) 
to describe the b quark in the decaying meson soft-collinear effective theory 

(SCET) to account for the energetic s quark in the final state 0, @, @, Gj], 11, 12 ]_, and non- 



relativistic QCD (NRQCD) to describe the quark pair in the charmonium final state 13l|-[l5||. The 
main adaptation of these theories to our case involves re-formulating NRQCD covariantly in an 
arbitrary frame (normally, NRQCD is formulated in the rest-frame of the charmonium state). The 
resulting effective theory, which we call covariant NRQCD (CNRQCD) accounts for the charm and 
anti-charm quark fields through four component spinors, to preserve covariance. 

We aim to show that the amplitudes for B decay to charmonium factorize. What is meant by 
this is that the amplitude for B — > ip + h, where ip is a charmonium state and h the light hadron, 
is the product of the amplitudes for B — > h + J\ and J2 — ► ip, where J\ : 2 are quark currents. That 
is, the B — ► ip + h is given by the form factor of the current J± between B and h states times 
the decay constant of the ip state. There is a simple physical picture which suggests factorization. 
Violations to factorization arise only if gluons are exchanged between the tp state and the B or h 
states. The effect of very energetic (hard) gluons affects the process at very short distances only, 
but will not affect the dynamical picture, so their effect can be absorbed into an overall constant 
coefficient. Moreover, this constant is calculable since hard gluons are perturbative. The effect 
of low energy gluons, on the other hand, is non-perturbative. However, the ip is a very compact 
bound state. In a multi-pole expansion, a long wavelength gluon interacts with the ip color-charge 
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distribution through it's color-dipole moment since the state is itself color neutral. But the dipole 
is very small because the ip is small. In the theoretical limit of very heavy charm, this coupling to 
the dipole vanishes. 

It should be noted that the weak interaction can also produce a charm-anti-charm pair with 
non-zero total color charge, in the octet configuration. The physical argument above indicates 
that in this configuration the quark pair does interact with soft gluons even at leading order. And 
while a color octet cannot form a (color neutral) charmonium state, the emission of a soft gluon 
can in principle turn the pair into a color neutral state. It would then seem that this electroweak 
contribution to the amplitude for ip production in B decays does not factorize. However, a long 
wavelength gluon interacting with a color octet leave the state in this octet configuration, and hence 
in a state which cannot produce a physical charmonium state. As soon as the gluon wavelength is 
short enough to discern the quark-anti-quark nature of the octet, a transition to the color singlet 
state can occur. But this is suppressed by the perturbative coupling constant of short wavelength 
gluons, so factorization holds at leading order and is expected to break down at the first order in 
the small expansion parameters. 

While this physical argument is not a proof, it does capture the basic ingredients needed to 
construct a more formal argument. As we will see, it will be important to separate the gluon 
degrees of freedom according to the wavelengths and frequencies with which different particles 
interact. This is accomplished by using the effective theory approach, combining HQET, SCET 
and CNRQCD. Degrees of freedom that are left out in this classification produce what amount 
to short distance corrections, which can be absorbed into re-definitions of coupling coefficients. 
On the other hand, active degrees of freedom exhibit some simplifications, as for example spin 
symmetries, that are of paramount importance in the solution to the problem. 

There is a quite extensive literature on the theory of B decays into charmonium. The naive 
factorization hypothesis in exclusive B — > ip + h, with tp and h denoting charmonium and strange 
states, respectively, has been tested against data, using models of the B to K form factor; see, e.g., 



Refs. lq. 1 1 71. ll8L 1 191] . Several theoretical approaches, distinct from ours, to factorization in QCD 



have been applied to exclusive B decays to charmonium. These include "perturbative QCD" [20j 



(pQCD) and "QCD factorization" [2ll . |22j. However, these approaches are not without trouble. 



Ref. [23(| shows that in QCD factorization infrared divergences arising from nonfactorizable vertex 
corrections invalidate factorization. And the pQCD approach does not find factorization but rather 
computes the non-factorizable part. 

One may also attempt to formulate factorization for inclusive charmonium production in B de- 
cays, as first proposed by Bjorken [24| (for semi-inclusive B — > J/tpX see 25]). This program was 
formulated more precisely in what is now also commonly known as NRQCD (13. [2a*|. This formu- 
lation, however, does not address factorization of amplitudes but rather of decay rates. Moreover, 
the program has run into difficulties as predictions of polarization of charmonium [27( have missed 
the experimental mark [2^ ]. Namely, the mechanism which resolves the factor of 30 discrepancy 
between theoretical predictions and experimental measurements of prompt ijj' production at the 
Tevatron involves gluon fragmentation to a sub-leading Fock component in the ip' wave-function 
and yields 100% transversely aligned ip n s to lowest order, which is not observed. The effective field 
theory approach presented in this paper may be used to study inclusive charmonium production 
in B decays. 

The paper is organized as follows. In section[H]we discuss the kinematics of the decay B — ► ip + h 
and outline the effective theory used to derive the decay Lagrangian. In section IIIII we list the 
degrees of freedom of the effective theory relevant in the leading logarithm approximation and give 
the Feynman rules necessary to do one-loop calculations in the effective theory. Also there we 
discuss the covariant modification of NRQCD. The effective theory electroweak Lagrangian in the 
leading logarithm approximation and at the leading order in the effective theory power expansion 
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is derived in section IIV1 In section [V] the simple formal proof of the factorization of the decay 
amplitude in the leading logarithm approximation is presented. In section IVT1 we perform the tree- 
level matching between the effective theory and full theory amplitudes with the soft (C)NRQCD 
gluons included and discuss the possibility of factorization beyond the leading logarithms. In 
Appendix [A] the differential decay rate for the semi-inclusive decay B — ► J/ip + h is evaluated. 
The initial values of the Wilson coefficients in the next-to-leading logarithm approximation at the 
matching scale of the decay are given in Appendix iBl 

II. KINEMATICS 

At the quark level the decay B — > ip + h can be described as the 6-quark decaying due to a 
weak interaction VF-boson exchange into a cc-pair and the s-quark. The spectator quark in the 
i?-meson becomes the second quark in h, see Fig. (pQ). 




u or d 



FIG. 1: The schematic representation of the decay. The heavy 6-quark, the heavy quark/antiquark cc-pair, 
the light s-quark, and the spectator quark are shown. 

The starting point of our analysis is the effective Lagrangian Cw of AB = 1 i?-decays of 
electroweak interactions obtained by decoupling heavy degrees of freedom of the SM such as the 
FF-, Z-boson and the i-quark at the scale \xw ~ M\y of the order of the FF-boson mass. Throughout 
we will refer to Cw as the full theory. The relevant part which governs at the quark level the decay 
b — ► s + (cc) is given by 

c w = -^v cb v: s (cMo + cMOz). (i) 

It has been obtained at the first order in perturbation theory in the electroweak coupling Gf 
(Fermi coupling) and the CKM-matrix elements V cb V* s . The flavor changing AB = 1 operators Oi 
are products of two quark currents with the spinor content [s . . . b] [c . . . c] which is the preferable 
choice when considering the formation of a cc-bound state as a final state. They read 

O = [sypP L b] [crfP L c] , 8 = [s lfl P L T a b] [c^P L T a c] . (2) 

These operators are related to the electroweak VF-boson exchange operators with the "natural" 
spinor content [s . . . c] [c . . . b] by a Fierz transformation. In the following we will frequently refer to 
Co and 08 as to the singlet and octet operators, respectively. 

The Wilson coefficients C{ depend on the scale /Uj ~ which is of the order of the 6-quark 
mass and contain the resummed leading QCD logarithms of the form ln(//b / fJ-w) to all orders in 
a s . Currently they are known up to the NNLO in QCD [2^ in a different operator basis though. 
The details of the transformation between two operator bases can be found in [2^] and the exact 
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definition of the evanescent operators involved in this transformation is given in Appendix [Cj Here 
we give only the numerical values of the coefficients 

C7^(M) = ^ + ^GfV)- (3) 

We find for the basis © at the NLO the central values C^ LO = 0.209 and C^ LO = 2.230 
with the input values 0: ot s (M z ) = 0.1176, M w = 80.403 GeV, M z = 91.1876 GeV, m b (m b ) = 
4.2 GeV (MS), /j,w = M\y and (i b = m b . It should be noted firstly, that the octet Wilson coefficient 
is one order larger than that one for the singlet and secondly, that the NLO correction to the singlet 
Wilson coefficient is roughly +50% since at LO: C^ 0) = 0.110 (whereas c| 0) = 2.228). Furthermore, 
the scale uncertainties amount to intervals C^ LO £ [0.137,0.269] and C^ LO £ [2.152,2.343] when 
varying the matching scale of electroweak interactions nw £ [55, 115] GeV and the scale fi b £ 
[2.5, 6.5] GeV. This variation is mainly due to fi b dependence which will be cancelled in physical 
observables once the matrix elements of the operators are calculated (up to a higher order a s 
residual scale dependence). 

Finally, at the scale [i b ~ m b the decay amplitude is proportional to the matrix element 

A(B ^ + h) ~ {tp, h\ C w \B) . (4) 

In the following we will set up the kinematics of the underlying partonic decay 6 — > s + (cc) in the 
full theory which will enable the construction of an effective theory in order to further investigate 
the structure of this matrix element. 

The wavelengths of the quarks participating in the decay are reasonably less than the wavelength 
associated with the confinement scale A c ~ 350 MeV -1 . The 6-quark and the cc-pair are considered 
to be heavy since their masses are roughly one order of magnitude and 3 — 4 times larger than 
the confinement scale, respectively. The s-quark is very light but it is also very energetic in the 
restframe of the 6-quark, so its wavelength is small. The numerical values of the quark masses are 
about 

m b {m b ) ~ 4.2 GeV, m c {m c ) ~ 1.25 GeV, m s ~ 0.1 GeV. (5) 

The masses of 6 and c quarks quoted are obtained from continuum determinations in the MS 
scheme [13]. The mass of the s-quark is not relevant for our calculation. We will work in the 
6-quark restframe (6RF) and decompose the momenta of the quarks accordingly. The 6-quark 
momentum can be decomposed into (omitting Lorentz indices) 

p b = m b [3b + h (6) 

with the scaling of the four velocitiy (3 b ~ 1 and the residual momentum k b ~ A c . 

We assume the cc-pair to have a small relative velocity v = \v\ <C c in their center of mass 
frame (cc CMF) in order to form the cc-bound state ip. In this kinematical configuration they are 
non-relativistic with a momentum scaling ~ m c (l + v 2 + O (v 4 ) ,v). This scaling persists also to 
an arbitrary boosted frame 

p c = m c p c c +Pc + k c , pc = m c Pcc -Pc +h (7) 

where the four veloctiy of the center of mass of the cc-pair scales as (3 C c ~ 1, the perpendicular 
components as pjr ~ m c v and the residual momenta k C) c ~ m c v 2 . Throughout the symbol v will 
denote the relative velocity of the cc-pair in their cc CMF whereas four velocities are denoted by 
Pi with i = {6, cc} and /3f = 1. 



The relative velocity v of the quarks of the cc-pair can be found from the self-consistency 
condition 



15|: 



v = a s (m c v), m c = m c (m c ) 



a s {m c v) 
a s (m c (m c )) 



(8) 



where the first equation follows from the virial theorem applied to the ground state of the hydrogen- 
like atom: Eq = — Ma 2 /2 = —Mv 2 /2 with M = m c /2 is the reduced mass. The second of Eqs. (|8j) 
is the LO expression for the running quark mass in QCD. For 5 flavors /3o = 23/3 and Aqqd = 0.225 
GeV. Numerical solution of Eqs. (JSj) gives: v ~ 0.61, m ~ 1.42 GeV, so that mv ~ 0.87 GeV and 
mv 2 ~ 0.53 GeV. The relative velocity v is not small although 2m ~ 2.84 GeV is close to the 
charmonium mass and mv 2 is comparable to A c . 

Since the effective theory is an expansion in powers of v around v = 0, the process b — > 
s + (cc) at the leading order in u-expansion is essentially a two-body decay and we can estimate the 
energy /momentum of the outgoing s-quark to be E s ~ 1.4 GeV in the 6RF. Then the energetic 
s-quark momentum can be decomposed as follows 



Ps 



(n-p s 



n I . .n 
+ P S + {n-Ps): 



(9) 



n) with n being a unit vector (the components are 

2 



n" = and (n • n) =2. Introducing p s = (n • p s ) 

J_ +\ i__ in i-x \ \2\ 

S 5 Ps J Ps 



E s (l, A, A 2 ) with A <C 1 when p s 



2 r S v « /2 

Here the two vectors = (1, n) and = (1, - 
given in the 6 RF) are light-like vectors, i.e. n 2 
and p^" = (n ■ p s ) the momentum p s = (pj ,p~t ,pf) scales as p 
is chosen to point approximately in the same direction as n. 

So, the wavelengths of all quarks involved are reasonably large compared to A c except for the 
spectator quark. Before the decay the spectator quark belongs to the l?-meson and after the decay 
it must become the constituent quark of h. 

Temporarily neglecting the spectator quark we will develop an effective theory for the matrix 
element in the limit where mj, — > oo, m c — > oo, and E s — > oo. We keep the ratio r = mc/m^ 
finite, otherwise the above limits are taken independently. The latter means that to develop the 
effective theory for the decay we can simply combine the leading orders of HQET (see [J, [1], 
and {HI), NRQCD [3], and SCET 0. We want to describe the decay in the 6RF, so we have to 
develop a covariant formulation of NRQCD. 



III. DYNAMICS: HQET, SCET, AND CNRQCD 

In this section we introduce the quark and gluon degrees of freedom of the effective theory and 
discuss a covariant form of the NRQCD. The Feynman rules necessary for the one- loop calculation 
in the effective theory are shown in Fig. We do not discuss HQET and SCET sectors of the 
effective theory and refer the reader to the extensive literature. 



A. Field operators 

In short, we need six fields to reproduce the infrared (IR) limit of the full theory amplitude: 

1. The outgoing s-quark is described by the SCET spinor i n ,p{x) (see e.g. jlj, and (33|). 
The large component of the collinear momentum p = \{n ■ p)n + p± becomes the label n, p 
on the field and the residual x-dependence is ultrasoft. 

2. The SCET collinear gluon An iP (x) transferring the collinear momentum. The momentum 
components scale like those of the collinear quark. 
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3. The incoming 6-quark is described by the HQET spinor hgAx) [34f | . 

4. The outgoing c- and c-quarks are described by CNRQCD spinors C/3 C c,p±( x ) an d i)^ Pi (2:), 
respectively. See the next subsection for the discussion. 

5. The ultrasoft field Au S (x). This field transfers momentum ~ m c t> 2 ~ E S X 2 ~ A c . 



(j.,a 



a) 



b) 



(P,k) 



np 



2 nk n p + p ± + iO 



rx^a 

ig s T n 1 



u.a 



c) 



: ig T a [#+ + t^±- 1 1 t^ f 
s n p n p' n pn p' 2 



d) 



Pb-k + iO 




(Pi, k) 



iPr 



01'- 



2 2 

pVk + ^+iO -(3 cc -k + ^+iO 



2m,. 



(l, a 



h) 



Pi 



(l, a 



= ±ig s T a |3^ 



Pi-qi 



H, a 



v,b 



-i gnv& 



|iv u ab 



qi + io 



j) 

k 

^' a vAAA/VW v ' b 



k) 



(q,k) 



r + io 



V,b 



i g^vSab 



n q n k + q^+iO 



FIG. 2: a) Propagator of collinear quark; n is the lightlike vector of the s-quark. b) Interaction vertex of 
the ultrasoft gluon and collinear quark, c) Interaction vertex of the collinear gluon and collinear quark, d) 
Propagator of the heavy 6-quark. e) Interaction vertex of the ultrasoft gluon and the heavy &-quark. f) 
Propagator of c-quark and/or c-antiquark in covariant NRQCD; f3 c5 is the 4-velocity of COM of cc pair. 
The second expression is the propagator of the charge conjugated field r) . g) & h) Interaction vertex of 
ultrasoft and/or potential gluon and c-quark in covariant NRQCD. Negative sign corresponds to the c-quark. 
i) Propagator of the potential gluon in Fcynman gauge. Momentum q± ~ m c v is transverse with respect 
to the 4-velocity of the cc-pair, /3 cE ■ q± = 0. j) & k) Propagators of the ultrasoft and collinear gluons in 
Feynman gauge. 
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SCET and NRQCD also include the soft gluons which transfer momenta ~ E S A and ~ m c v, 
respectively. At one loop the soft gluons do not contribute to the effective theory decay amplitude, 
so they're not listed above. The soft gluons are discussed in section IVII where tree- level matching 
to all orders in g s is performed. 



B. CNRQCD 

1. Lagrangian 

In the center-of-mass frame of the cc pair (cc CMF) the effective Lagrangian is provided by 
NRQCD. The desired covariant Lagrangian can be obtained from the NRQCD Lagrangian [35[ 
in the cc CMF by a Lorentz boost to the frame where the center-of-mass velocity is p cc -. At the 
leading order in u-expansion this gives: 



Pi 



Pec ■ iD + 



(pi 



2m c 



t- —C 



PA 



Ml! 

2m c 



Pec ■ iD 



Cpx- (10) 



Here the covariant derivative D contains only the ultrasoft field A us . When calculating one-loop 
QCD corrections to the effective theory Lagrangian it is more convenient to work with the charge 
conjugated field rf^ c _ p± because it creates the anti-quark c in the final state and this is exactly 
what we need in the effective theory for the B — > ip + h decay. 

At the leading order in ^-expansion the effective theory Dirac spinors £/3 c5iP± and rj^_ p ^ are 
related to the full theory field operator ip as: 

ij,{x) = e-im^-x J- e-^*Zp^ Px {x) + e im ^ x e^^Jx) + O (v) , (11) 

with p± the covariant potential labels. We use symbols £ and ij for the quark and anti-quark 
fields following the conventional NRQCD notation. However unlike in the standard approach we 
keep all four components of the spinors to ensure relativistic covariance. The 4-spinors satisfy the 
constraints 

p pJp^ = °> Ppjife^ = °- ( 12 ) 

Here = (l+@)/2. The solutions of (|12p are such that in the cc CMF their upper components be- 
come the conventional NRQCD two-component Pauli spinors £ p and rj p and the lower components 
vanish: 





f/Wx = V I ftffl c ' VPcc,p± ~ \ 7T^ I p^ff „ • ( 13 ) 



In other words £/3 cS .p ± and rj/3 c5tP± are the NRQCD spinors £ p and rj p boosted into the frame 
with velocity — /3 C c(/3c 5 ) _1 with respect to the cc CMF. In the cc CMF p c5 = (1, 0) and the standard 
form of the NRQCD Lagrangian is reproduced. This formalism is essentially the same as the one 
discussed in [if]]. The spinor rf^ a _ p± satisfies 

P f3 c ^Zc, P± = VPce, P± ' ( 14 ) 

where Pa = (1 — $)/2, and transforms under the same representation of SU(3) as £p cS>p± - The 



explicit form of ffp c - p± in terms of Pauli spinor rj p is given by ([T 



9 



NRQCD also includes the off-shell quark modes called soft quarks. The energy and momentum 
of the soft quark both scale like ~ m c v and therefore the Lagrangian for the soft quark at the 
leading order in ^-expansion is the same as HQET Lagrangian. However we will see in section IIVI 
that at one loop the soft gluons do not contribute to the decay amplitude and so we do not discuss 
them here. 

The NRQCD part of the decay amplitude also includes a 1/v (Coulomb) singularity. This 
singularity is reproduced by the local (on the ultrasoft scale) four-quark operator (see [33] and [151 ]) 
corresponding to the Coulomb interaction between the heavy quarks. In this paper instead we 
have adopted a simplistic approach: we introduce the potential gluon that transfers the frequency 
~ m c v 2 and momentum ~ m c v (in the cc CMF) and does not remove the quarks off-shell 35f | . 
This is enough to reproduce the Coulomb singularity at one loop. The Lagrangian that couples 
the potential gluon and the on-shell quarks changes the momentum label p± on the quark fields £ 
and rf and gives rise to the same interaction vertex as the ultrasoft gluon. 



2. Full theory spinors in terms of £ and rj 



Matching a full QCD matrix element onto the matrix element of the CNRQCD requires writing 
down the full QCD spinors u and v in terms of the CNRQCD spinors £ and rj. This is done as in 
NRQCD only the frame of reference is not specified. For example: 



m)u 



-2mP g + $ c )u 



0. 
0. 



(_ 2m p-+p- ? j c )(P+ + P / 



P 

-2mP; 



P 



P^u 



P 



P 

P-PcP, 



U- 



+ P pi ) c± 

(2m + (3 ■ p c )u- + Pp$c±€p, 
-Pptfc±£p c5 ,p ± , 



P 



u 



PX 



2m + j3 ■ p c 
1 

1 + 



2m + (3 ■ p, 



£.Pc 



= o, 

U = CPcc,PJ 

= 0, 



/3=(1,0) 



■PA 



JP 



(15) 



E+m^P 

Here p c denotes the sum of the perpendicular and residual components in the momentum decompo- 
sition given in 0. In the equation the effective theory spinor £/3 cfiPx is introduced as a projection 
of the full QCD spinor u onto the subspace restricted by the constraint (|12p . The second half of 
the bispinor subspace u_ accounts for the difference of the order v between the full and effective 
theories. One can see that in the cc CMF the full QCD spinor is reproduced. Unlike the £-spinor 
which has only two upper components in the cc CMF, the full QCD spinor has also two non-zero 
lower components suppressed by the first power of velocity in the cc CMF. 

To parameterize the v-spinor of the full theory by the rj '-'-spinor notice that the charge conju- 
gated spinor v c = Cv* satisfies the same equation as the ?x-spinor. The CNRQCD Lagrangian for 
the anti-quark is derived from the full theory using charge conjugation. Therefore the t> c -spinor 
must be written in terms of the effective theory 77-spinor in the same way as the w-spinor is written 
in terms of the £-spinor and then the v-spinor follows after charge conjugation: 

f 1 1* 

^ 1 + „- I a _ P/3 fic± rffa, 



■ 2 c* 

27 V 



1 



2m + (3 ■ p, 

Ptic± 



PX ' 



2m + (3 ■ p c 
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1 



2m + (3 ■ p c 



Cpx" ( 16 ) 



Here, componentwise: 



Cpx = (^Vh, P J = ( fjp ] where (-ioV) = (-^{f. (17) 

Notice that the spinor rj^_ p ^ satisfies (|14p . i.e. for the charge conjugated spinor the 4- velocity is 
reversed. 



IV. EWET LAGRANGIAN IN THE LEADING LOGARITHM APPROXIMATION 

To obtain the effective theory electroweak (EWET) Lagrangian we perform the tree-level match- 
ing between the full theory amplitude of b — > (cc)s decay and the decay amplitude in the effective 
theory to all orders in g s and at the leading order in v and A. Then the result is RG-improved 
using the anomalous dimension matrix extracted from the results of one-loop matching which cor- 
responds to summing the leading Sudakov logarithms and gives the EWET Lagrangian in the 
leading logarithm approximation. The soft gluons do not contribute at one loop and therefore do 
not modify the anomalous dimension matrix in the leading logarithm approximation. From this 
we infer that in this approximation the soft gluons can be ignored. 

A. Tree-level matching 

At tree level the matching is done by replacing the full theory fields of ([1]) with the effective 
theory fields of section IIII Al and using the reduction formula (see Appendix [D]) 

I^Pl ® I^Pl -» \pr ®4±~ \(Pcc • n)P R <g> 75 - ® U± (18) 

to decompose the full theory Dirac structure into the sum of the three Dirac structures of the 
effective theory. As we have inferred the collinear gluons are the only degrees of freedom giving 
rise to the local operators on the effective theory side in the leading logarithm approximation. 
Then the result of tree-level matching follows easily from the O (g s ) calculation and the SCET 
invariance of the effective theory Lagrangian that requires the collinear gluons and the collinear 
quark fields to combine in the jet field [||: 

G (tree) = [inj.WTiCh^fo^TjC^.pJ. (19) 

Here C = 1 or T a and Tj (8) Tj = {Pa <8> Pr ® 75, ^- <8> The sum over potential label p± 

is understood. The operator W is the Wilson line 

w = E ex p i~3'^p n ■ A ^i) ( 2 °) 



perm 



where V is the operator that picks up the net label momentum on the right. Momentum con- 
servation implies that at each order in g s the labels on the collinear fields satisfy p + Y2i Qi = 
rribPb — 2m c (3 C c- The corresponding Feynman diagram in the effective theory is shown in Fig. [3j 
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FIG. 3: The effective theory amplitude at the tree-level. Double line corresponds to the incoming heavy 
6-quark (labeled by 4- velocity 0b), the dashed line corresponds to the outgoing light s-quark (labeled by 
collinear momentum p) , the solid lines correspond to the outgoing quarks of cc-pair (labeled by the potential 
momenta ±p±). For the c anti-quark the fermion flow is reversed as in the full theory. Only one collinear 
gluon incoming to the electro- weak vertex is shown. 



B. Decay amplitude at O (a s ) 

Here we present the results of one- loop calculation of the decay amplitude in the effective and full 
theories. Feynman diagrams are calculated in the MS-scheme. The momenta of the four external 
particles are set on-shell and the corresponding IR divergences are regulated by the dimensional 
regularization. 

1. Effective theory 

Overall in the effective theory there are eight amputated diagrams. Seven of them are due to 
loop corrections to the zero order term in the g s expansion of the tree-level Lagrangian (|19f) : 

Off = [in^CihpMp^TjCivg^], (21) 

where i = 0, 8 stands for the singlet and octet operators and j = 1,2,3 for the Dirac structures 
introduced in (fT8j) . 

Six diagrams where the ultrasoft gluon is exchanged between the quarks are shown in Fig. [H 
Two more diagrams in Fig. [5] correspond to the potential gluon exchange between the quarks of 




FIG. 4: Diagrams with the ultrasoft gluon exchange in the effective theory. 

the (cc)-pair and the emission of collinear gluon from the electro-weak vertex. 

The diagrams involving the soft CNRQCD gluons emitted from the electro-weak vertex are 
shown in Fig. [6j All these diagrams include also the off-shell soft quarks. The propagators for the 
soft quarks are HQET propagators because the term /3 c5 • q s ~ mv and the kinetic energy term 
in the CNRQCD propagator p^_/2m ~ mv 2 must be discarded. This makes Feynman integrals 
corresponding to the diagrams in Fig. [6] scaleless and thereby identically zero in the dimensional 
regularization. So the soft gluons do not contribute to the amplitude at one loop. 
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FIG. 5: The potential gluon exchange between the quarks of cc-pair and collincar gluon emission from the 
electro- weak vertex. 






FIG. 6: Soft gluon exchange between the quarks of cc-pair and soft gluon emission from the electro-weak 
vertex. The zigzag line corresponds to soft gluon. The double lines for c and c-quarks stand for the soft 
propagators. These diagrams are identically zero in the dimensional regularization. 



The non-renormalized amplitude is given by the sum of the diagrams in Figs. U] and The 
ultraviolet (UV) divergences are then removed by the field renormalization factors for the heavy 
and light quarks, Zu and Zi, respectively, and by the renormalization matrix of the electroweak 
operators for i,j = 0, 8, Z\j\ 
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In 
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+ 
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fl 



In 



2 euv 
6 



+ 
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euv 2E S 



Zos - 
2m c 
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Z 



80 



In 



+ 



3vri 19 



Suv 



+ 



0. 
1 



3 euv 



(22) 



Here E* 



-Ami 



is the energy of the s-quark in the 6RF. Finally the MS-renormalized amplitude 
is multiplied by the Lehmann-Symanzik-Zimmermann (LSZ)-factor \fR. The latter is computed 
from the residues at the poles of heavy and light quark propagators also renormalized in MS-scheme. 
This gives the on-shell matrix element: 



The resulting expression for the amplitude is given by: 

i=0,S i=l,2,3 



where 



(23) 



(24) 
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a s 

47T 


[C F 
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(25) 
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At one loop there is no mixing between the singlet and the octet operators: Zq$ = Z^o = 0. One 
can see why by writing down the color structure of the effective theory diagrams. The diagrams 
of the first column in Fig. H] and the diagrams in Fig. [5] do not alter the color structure of the 
vertex. The remaining four diagrams in Fig. |4] actually mix the singlet and octet operators but the 
interaction vertices for the c and c quarks have different signs as it follows from the Feynman rules 
of the effective theory and therefore the color structure of the sum of the second and the third 
and /or the fifth and the sixth diagrams in Fig. 0] is either 

rC®[T a ,C] or CT a ®[T a ,C}. (26) 

For C = 1 the commutator vanishes and for C = T a the direct product is reduced to T a ®T a . 



2. Full theory 

In the full theory the nonrenormalized amplitude is given by the sum of six diagrams in Fig. 
We are using the NDR scheme for 75. The on-shell value of the full theory amplitude follows when 




FIG. 7: Diagrams contributing to the amplitude b — > (cc) s in the full theory. 

the sum of the six diagrams is multiplied by the field renormalization factors Z q , the operator 
renormalization matrix Cf s,Te = CjZji, and the LSZ factor: 

7-1 a °r 1 7-s^ asl ( 06 \ 

— 1 ~ 1~°F , ^ij — Oij -+- 4 , 

4tt euy 4vre uv \ 3 ~ 2 / 

V^^s = l + ^f-^--21n^4-41n^-8^ . (27) 
47r V 3euv ml m l c ) 

The final step is projecting the full theory Dirac structures to the Dirac structures of the effective 
theory which corresponds to taking the limit v — > and A — ► 0. This is done by means of the 
reduction formulae listed in Appendix [Dj The projected full theory amplitude contains only the 
Dirac structures of Eq. (fl~8j) . 



3. Matching 

At one loop the matrix element on the full theory side upon expanding the Dirac structures in 
terms of the effective theory structures becomes: 

[Af + ^m). (28) 
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Here the index i is cumulative and labels both the Dirac and color structure. On the effective 
theory side the matrix element has the form: 



C (0) 0^(1) 

1 4tt 1 



<<*> + £*g } <°i) 



(29) 



where C- ^ is the first order correction to the Wilson coefficient and the second bracket is the matrix 
element in the effective theory at O (a s ) as it comes from the one-loop diagrams. Equating (j28H 
and (f29j) gives the matching condition at one- loop: 



C 



The coefficients cf^ have been introduced in the tree-level matching formula (j 18() : 



(o) 



.4 



(o) 



and C} 



(i) 



.4 



(i) 



.4 



(i) 



AO) y(l) 



(30) 



Cf\r) G 



1 1 

\~r '~2 



with 



(31) 



The analytic expressions of the coefficients C\ are given in Appendix [Bj Their numerical values 
at the matching scale of the decay are discussed in section IVI Bl 

The coefficients are given by Eq. (I25|) and contain the IR divergences of the effective theory. 
The coefficients contain both the finite part and the IR divergences of the full theory. We 
confirm the cancellation of the IR-divergences between the full and effective theory amplitudes. 
From this fact we infer that at the leading order in the effective theory power expansion and in 
the leading logarithm approximation the effective degrees of freedom include only those listed in 
section lill Al This is the main result of the paper. 



C. Summing the leading logarithms 

The renormalization matrix Zy in Eq. (I22p is diagonal and contains the terms ~ 1/e 2 , ~ 
(1/e) log /z, and ~ 1/e. The terms proportional to 1/e 2 and ~ (1/e) log /i give rise to the leading 
Sudakov logarithms ~ a" ln n+1 (u) as we will see shortly. The terms ~ 1/e give rise to the next- 
to-leading logarithms ~ a™ ln n (^) and should be combined with the next-to- leading logarithms 
coming from ~ 1/e 2 terms in the renormalization matrix Zij calculated at two loops (see Table II 
in 0] for the counting scheme of the divergent terms). 

Discarding the terms ~ 1/e gives the renormalization factors Zqq = Z^s which turn out to be 
equal to the renormalization factor for the (s6)-current (see 0]): 
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euv 2£/ s 



(32) 



Therefore the anomalous dimension matrix (ADM) in the leading logarithm approximation is the 
same for all operators Oij and coincides with the ADM of the (sfr)-current in SCET whose Wilson 
coefficient is calculated in 0j: 



C 



LL 



2E S 



C 7 (0) (r)exp 



oo 

E 

n=l 



^£_F ^ ("I) 
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2E X 



(33) 



We have chosen the form (|33|) for the Wilson coefficient to make it explicit that solving the RG- 
equation is indeed equivalent to summing up the IR double Sudakov logarithms. 
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SCET gauge invariance of the jet field [£ njP W] ensures that loop corrections to this operator and 
to £n,p are the same 0] . This observation allows one to combine the tree- level matching formula (|19h 
with the Wilson coefficient ([33]) by treating the latter as the operator function of n • V* inserted 
on the right to the jet field [9]: 

V i=0,8 j=l ^ ^ ' 

(34) 

Here picks up the net label momentum of the collinear operators on the left and Ci(//&) are 
Wilson coefficients ([3]). The operator n ■ V' picks up the total momentum of the jet at each order 
in the g s expansion. At tree-level this amounts to replacing 2E S — > n ■ but when evaluating 
a s -corrections the presence of the operator function in (|34p becomes nontrivial. The dependence 
of C\y on the scales [i and fi b is to remind that the Lagrangian is obtained in two steps: firstly by 
running from the electroweak scale to /z& = 2E S and then to fx. 



V. FACTORIZATION 

Applying the EWET Lagrangian (|34h to the decay of the l?-meson requires evaluating the 
matrix elements of the Lagrangian between the hadronic states of the effective theory. In this 
section we prove that in the leading logarithm approximation the matrix element of the EWET 
Lagrangian between the effective theory states factorizes into the product of the matrix elements 
of the heavy-to-light current and decay constant of the charmonium state. 



A. Ultrasoft-free EWET Lagrangian 

As we have seen the ultrasoft gluons are the only degrees of freedom common to all three sectors 
of the effective theory in the leading logarithm approximation. The following field redefinition (see 
e.g. makes it possible to rewrite the EWET Lagrangian (|34l) in terms of the ultrasoft-free field 
operators by incorporating the ultrasoft gluons into Wilson lines: 

c — y t(0) a _y /(0)yt w — y w(0) yt 

ht = Y p b h f b , 

c v AO) C v 0(0) / OK x 

Sf3 ceP± = Y /3ce£f3 c5 ,p ± > V(3 C c-p ± ~ Y PocVf3 ce -p ± - (30) 

Here Yj, with I = n, fib, or /3 c5 , stands for the path ordered exponent: 

Yi(x) = Pexp(ig, J dsl ■ A us (l s)J . (36) 

It is a well-known fact that ultrasoft gluons do not renormalize the four-quark Coulomb interaction 
term in NRQCD (see e.g. [3]) and so, as long as soft gluons are ignored, the ultrasoft-free CNRQCD 
Lagrangian includes only the ultrasoft-free potential quarks and potential gluons. 
Applying these field redefinitions to (f34"|) gives: 

Of = &W®C? r) Y^ jCi Y Pb ^[^^LriW^^J- ( 37 ) 

It is instructive to pause and check that at one loop corrections to (|37|) exactly reproduce the 
effective theory amplitude given by ([23]) and (f25|) . Expanding the ultrasoft and collinear Wilson 
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lines in (|37p to the first order in g s reproduces the diagrams shown in Fig. H] and the second diagram 
in Fig. Terms of the order gl coming from expanding the Yi(x) are self-energy diagrams and 
must be discarded since only the amputated diagrams contribute according to LSZ prescription. 
Wilson lines Y[{x) bring the renormalization factors which reproduce the renormalization factors 
of the quark fields due to the ultrasoft self-energy diagrams. Expanding Yp c - reproduces only the 
HQET part of the potential quark propagator in Fig. [21 This is sufficient because the corresponding 
diagrams (2nd, 3rd, 5th, and 6th diagrams in Fig. [3]) do not have 1/u-divergences. Physically for 
the ultrasoft gluon the CNRQCD quark is the same as HQET quark as long as the relative motion 
of the quarks in the cc CMF is ignored. 

B. Matrix element for B — > ipX s . 

Consider the matrix element of (f37l) between the states (X sn , ip/3 cE \ and \0/3 c5 , Bp b ), where X sn 
is any state which includes a collinear s-quark. The states |0^ c -) and {ipp c7i \ do not contain ultrasoft 
gluons, so the matrix element of the singlet operator factorizes: 

(X sn \ &°> (0) Cf L r) r,y n t Y Pb hf h \ \B Pb ) (^J [$l p± T 3 v l%J |0 fe ) . (38) 

Now let us show that the matrix element of the octet operator (C = T a ) vanishes and therefore 
the color state of the cc-pair does not contribute to the decay rate (in the leading logarithm 
approximation). To this end it would suffice to show vanishing of the correlator 

(X sn , ^{0 { S{*)&W^Cf L r) T jY ^Y Pb hf b ] 

x fSLr, Ylj^ n^%J(y)} |0 fc , B Pb ) , (39) 

where 0^)(x) is the ultrasoft-free interpolating operator carrying quantum numbers of tpf3 ci ,, the 
bound state of the cc-pair. The matrix elements (|39|) and (|38|) with Cj = T a are related by LSZ 
reduction formula and vanishing (|39]) implies that (|38|) is also zero. 

The diagrams contributing to the T-product in (J39]) are the ladder diagrams shown in Fig. [HI 
where the potential gluons are exchanged between the quarks of cc-pair. The interpolating operator 




FIG. 8: Ladder diagrams that contribute to the L-product ([39]) at the leading order in a s . Wiggly lines 
represent ultrasoft gluons which couple between themselves and to the external state B. The collinear gluons 
(not shown) couple to the collinear s-quark line and to the collinear degrees of freedom in X s . 
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(x) and the currents in (|39j) are color singlets, therefore the color structure of the diagram in 
Fig. El is: 

^ j^l, ai,a2---'p r l^^ 1 T b2 [ [T a T ai ] T a2 ] ] T b ' 2 T bl | (40) 

Here o« and bi are adjoint indices. The commutators inside the trace arise upon expansion of the 
ultrasoft operators Yp c -. The insertions of T bi on both sides are due to potential gluon exchange. 
The quantity A a ' ai ' a2 represents the color structure arising when contracting the ultrasoft fields 
from the cc-current with the fields of the heavy-to-light current, and with the ultrasoft degrees of 
freedom in the Bp b state. Expanding the commutator inside the trace gives the structure propor- 
tional to T a . Contracting the indices bi on both sides of T a again gives a structure proportional 
to T a . Therefore the whole trace is proportional to the trace of T a and vanishes. 



VI. BEYOND THE LEADING LOGARITHMS 



If it were not for the soft NRQCD gluons next-to-leading corrections could be included sys- 
tematically into Cij(n ■ V* / V, r) of Lagrangian ([34"]) . Without the soft gluons we simply have 
three HQET quarks, two of them propagate along the same Wilson line in the opposite directions 
and as the matching shows decouple from the (sb) current. To see explicitly whether soft gluons 
violate factorization in the next-to-leading logarithm approximation two loop calculation in the 
effective theory is required. We have not attempted this calculation, although made the first step 
in this direction by performing the tree-level matching for the EWET Lagrangian that includes 
both collinear and soft NRQCD gluons. 



A. Including soft gluons 



The matching is performed using the technique explained in the Appendix A of [32]. The method 
is to introduce the auxiliary fields corresponding to the off-shell quarks and gluons, to write down 
the Lagrangian for the auxiliary and on-shell fields, and then to integrate out the auxiliary fields 
by solving the equations of motion (EOM). The integration corresponds to the tree-level matching 
to all orders in g s . The result given by (|49p does not exactly reproduce the corresponding result for 
the heavy-to-light current in SCET (see e.g. Eq. (A22) in [32|) because the soft NRQCD gluons 
are not protected by the SCET symmetry under the soft gauge transformations. In this section we 
briefly outline the derivation. 



1. Auxiliary fields 



The first step is to draw all tree-level diagrams that introduce the off-shell fields. Figures [9] 
and [10] show how the off-shell modes arise for the quarks of cc-pair and s-quark, respectively. 

For the 6-quark we have to match in two steps. Both soft and collinear gluons give the 6-quark 
an off-shell momentum but the corresponding off-shellnesses are of different orders of magnitude. 
A soft momentum is ~ mv and collinear momentum is ~ m. Once the collinear momentum 
is transferred to the 6-quark any off-shell momentum brought by the soft gluons is subleading 
compared to the collinear momentum propagating along the line. A diagram where a soft gluon is 
inserted between collinear gluon(s) and the electro-weak vertex contributes the factor 
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FIG. 9: a) Collinear gluon interacting with c-quark and the off-shell mode c'. Interaction between collinear 
gluon and the off-shell mode c' is given by the same vertex, b) Collinear gluon interacting with ?7 c -field of 
the antiquark and the off-shell mode rj c/ . Interaction between collinear gluon and the off-shell mode rj c/ is 
given by the same vertex. The overall sign is opposite relative to a) due to the reversed fermion flow. 
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nq s 



pf-q s vt 



FIG. 10: Soft gluon interacting with s-quark and the off-shell s'-ficld. Interaction between collinear gluon 
and the off-shell mode s' is given by the same vertex. 



where q c is the collinear momentum propagating along the 6-quark line. (A gauge field scales 
like the momentum it transfers.) Therefore only the diagrams with the soft gluons next to the 
heavy quark field hp b and the collinear gluons next to the electro- weak vertex shown in Fig. 1111 will 
contribute. To account for this we introduce two off-shell fields for the 6-quark: b' which is off-shell 




FIG. 11: Diagrams which contribute to the tree-level matching of the heavy 6-quark. Soft gluons (zigzag 
lines) must be next to the heavy quark field. 

by ~ mv and b" which is off-shell by ~ m, (see Fig. [T2j) . 

In addition to the soft and collinear gluons we should take into account the mode arising when a 
soft gluon fuses with a collinear one, see Fig. [T3l The auxiliary gluon field Ax transfers momentum 
qx ~ m(l, v, v) in the light-cone coordinates (compare to the collinear momentum q c ~ (1, A, A 2 )). 
Recall that v ~ A for the decay kinematics. At the leading order in the effective theory power 
expansion the interaction vertex of the field Ax with the heavy quarks is the same as for collinear 
gluons. The interaction vertex of Ax and the s-quark is the same as for a soft gluon. The gluon 
Lagrangian for Ax is discussed in [32]. 



2. Matching 



Writing down the auxiliary Lagrangians which generate the Feynman rules listed above is 
straightforward. The EWET operator is given by: 

(trcc) = \(£n, P + s)T 3 C (b" + b')][(^ cP± + c')T,C {if + (42) 
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FIG. 12: a) Soft gluon interacting with the 6-quark and the off-shell mode b 1 . The interaction between the 
soft gluon and the off-shell mode b' is given by the same vertex, b) Collincar gluon interacting with the 
off-shell mode b' and the off-shell mode b". The interaction between the collinear gluon and the off-shell 
mode b" is given by the same vertex. 




FIG. 13: Soft gluon fusing with a collinear gluon into the off-shell field Ax- 



Here the auxiliary fields are understood as perturbative expansions of solutions for the correspond- 
ing EOM following from the auxiliary Lagrangians. For the 6-quark we have to solve first the EOM 
for the b' field and then use the solution to solve EOM for b". The solutions are 

c =Cwx04-i) ; 

s' = £n,p(S nX — 1), 

b" = {W x - l)b', 

b' = Sp b h Pb , (43) 
where Sp b , Wx, and Sx are Wilson lines 

s p b = Yl exp \7r~p 9sl3b ' As >i) ' 

perm ^ ' 

W x = exp (~-^p9s[n ■ A x , q + n- A C)9 ] J , 

perm ^ ' 

S lx = exp \9s[n-Ax, q + n-A Stq }-^—] . 



perm 



(44) 



Combining Eqs. ([52]) . fljgj) , and (jUJ) gives: 

(tree) = [ZnfSlcWxTjC SpMfa^TiCfa^J. (45) 

Solving the EOM for the auxiliary field Ax makes it possible to write 0( tree ) in terms of A c 
and A s only. This remarkable calculation is performed in Appendix A of [32|], to which the reader 
is referred. The only difference is that here we work with the soft NRQCD gluons instead of the 
soft SCET gluons. The result is: 

S x W x = WSl (46) 
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where W is the collinear Wilson line (|20p and S n is given by: 

perm ^ ' 

Using (|46l) eliminates the auxiliary field Ax from (|45p which gives: 

0(**o = [e niP w4riC5 A /i A ][e/9 ce p x r i c^ ce _ p j. (48) 

It is easy to verify this formula at O (g s ) by explicit matching. 

Now the net momentum of the cc-pair is transversal with respect to the 4- velocity f3 c5 . Therefore 
the net momentum q s of the soft gluons emitted from the EWET vertex must be transversal as 
well because they end up inside the (cc)-current: (5 C c ■ q s = 0. This observation suggests that (|48j) 
must be written as 

o (trcc) = [Zn*W6te. v & n r j cShhh][tfa Px r j c v g B _ p ^ (49) 

where the Kronecker delta enforces transversality of the net soft momentum transferred from sb- 
to cc-current. Equation (|49p is the final result. 



3. A remark 

Equation (|49p does not hint that the factorization holds in the next-to-leading logarithms ap- 
proximation and beyond. However, the transversality constraint results in something resembling 
the factorization. Consider Coulomb gauge. In this gauge gluons are transversal, so the soft gluons 
obey: /3 c g • A s = 0. Now let us apply the transversality constraint in (|49p separately to the soft 
momenta emitted from the b- and s-quark lines. (Actually the constraint restricts the net momen- 
tum.) Then the Wilson lines S n and Sp h become equal. To see this notice that at the leading order 
in power expansion the momentum conservation reads: 

2m c E s 

Pb = Pec H n. (50) 

Applying this equation to the argument of Sg b (|44p together with the transversality constraints 
yields: 



Pb ' A S) q Tl ■ A s< q 



b -V n-V 



Sn = Sp b . (51) 



The last equation allows one to move the soft gluons from the (sb) current to (cc)- current, so 
the factorization "holds" even when the soft NRQCD gluons are included. However the attempts 
to make anything rigorous out of the above observation were unsuccessful. Most probably because 
the factorization is broken beyond the leading logarithm approximation. 

B. Wilson coefficients in the next-to-leading logarithm approximation at \i = rrib 

The one-loop matching (see section IIV B 3p gives the initial values of the Wilson coefficients 
at the matching scale [i = fib- At this scale the EWET Lagrangian in the next-to-leading 
logarithm approximation is: 

C^ b ,, b ) = -^V cbV ; s £ Ea(^){c'f(r)4 fc + ^cfg,.(^ 1 r)}og.(52) 

v i,k=0,8j=l ^ b 
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The numerical values of C\^. - for ^ = nib an d r = m c /mb = 0.30 are: 

(i) / -11.27 -11.27 - 3.41A 

°ik;i^' - ^ -2.50 - 0.76i 8.94-11.55* )' 

r (i) n n <un- f 1611 19M + 6 - 58i ^ 

°ifc;2^' u -°^ ^ 4.36 + L46i _13.56 + 20.29i J ' 

(1 30) - ( " 10 - 75 10 - 81 + L23i ^ ( 5S ) 
^ik;s^ U-M) ~ [ 2.4o + .27i -7.87 + 9.00i J ' [b6) 

Here indices i and A; label the rows and columns, respectively. 

Comparatively large numerical values of the non-diagonal elements i ^ k of the coefficients 
^ik-j together with the observation that numerically the Wilson coefficient Cg LO (rrib) is about ten 
times larger than the Wilson coefficient of the singlet (see section UTJ) suggest that the next to the 
leading logarithm correction to (|34p is of the same order of magnitude. 



VII. CONCLUSION 

In this paper we have proposed the effective theory for the decay B — > ip + h, where tp is a 
charmonium state and h a light hadron neglecting the effects due to the spectator quark. We have 
identified the relevant degrees of freedom and verified explicitly by doing the matching calculation 
at one loop that the effective theory reproduces all the IR divergences of the perturbative QCD. We 
have also shown that in the leading logarithm approximation the effective theory decay amplitude 
factorizes into the decay constant of the charmonium state and the matrix element of heavy-to-light 
current. The latter has been extensively studied. 

However including into consideration the soft NRQCD gluons most probably ruins the fac- 
torization beyond the leading logarithms. At least the tree-level matching result (I49h does not 
suggest the opposite. The comparatively large value of a s ~ v ~ 0.6 in charmonium implies that 
the higher order corrections cannot be ignored. To achieve a better quantitative understanding 
of non-leptonic .B-decays into charmonium we should definitely go beyond the leading logarithms. 
Our results can be used e.g. for calculating the RG-improved Lagrangian in the next-to-leading 
logarithm approximation provided the anomalous dimension matrix of the effective theory is known 
to two loops. 
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and to prof. A. Manohar for discussing the results. 
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APPENDIX A: SEMI-INCLUSIVE DECAY B -> J/ip + h. 

In this section we apply the EWET Lagrangian (|34p to the calculation of the semi-inclusive 
decay rate of B — > J/ifj + h. Due to factorization of the charmonium state in the leading logarithm 
approximation the decay can be treated along the same lines as the semi-leptonic decay B — > 
[lOfj. The derivation below is presented in some detail. The result is given by (1A20|) . 
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1. Kinematics 



The differential decay rate is given by 

{J/4>, X S \H W \B)^ 2 



dr 



dE 



\Pj/4'\ \ - 



(2ir) 4 5 4 (p B -pj/ip -px)- 



2M B 



(Al) 



Here \pj/$\ 
range 



E 2 



M 2 



J/1> 



. The differential decay rate depends only on Ejj^ that varies in the 



E 



J/1> 



M 



M B + M 



J/1» 



J/V 



2Mb 



[3.1, 3.56] GeV. 



(A2) 



There is a region of Ej/^ where our effective theory can be used. In this region the state X s must 
be collinear, i.e. the SCET parameter A = Mx/2Ex, where Mx is the invariant mass of the jet, 
is small. On the other hand Mx must be reasonably larger than A c , so the jet is inclusive. 

It is customary to take A ~ ^K c /Ex <C 1 as a value of SCET parameter for a jet which is both 
collinear and inclusive. The kinematic expressions for Ex and Mx in terms of Ej/^ are 



Ex = Mb — E 



and 



so we have to look for a region in (1A2I) where 

M x _ 
2E X " 



Ml 



M 2 B + M^-2M B E J/i ,, 



(A3) 



Ex 



< 1. 



For A c ~ 0.35 GeV such a region exists near the center of the range 
Mx ~ 1.6 GeV. The SCET parameter is reasonably small, on a par with v 



(A4) 

2J where A ~ 0.4 and 
0.6 in the charmonium. 



2. J/ip decay constant 

The matrix element in (|Aip factorizes according to (|38|) . So, the first step is the evaluation of 
the matrix element of the (cc)-current (see (|18p ): 

(J/fyJ [4!Lx r ^fe- P x] l°A*> ' where r i = {^»75,^ + ±}. (A5) 

The matrix element of 75 vanishes by parity; it is non-zero for the pseudoscalar state r\ c - The other 
two can be expressed via the J ftp decay constant. In the full theory it is defined as (see (38|): 

<JMp, e)| mx)^(x)] |0) = -if m M m e^(p) e**. (A6) 

Here Mjj^ is the mass, and £*Ap) is the polarization vector of the state. The decay constant 
determined from the leptonic decay mode J/t/} — > l + l~ is fji^ = (405 ± 15)MeV 

To match the left-hand side of (|A6p to the CNRQCD at the leading order in u-expansion we 
replace the full theory operators with their expressions in terms of the ultrasoft-free effective theory 
fields (see (fTT|) and ([35]) ). The states of the full theory are replaced with the states of ultrasoft-free 
CNRQCD. In the sum over p± and p±' only the terms where p± +p±' = survive. Also we should 
take into account the Wilson coefficient Cj/^(m ret i/n), where m re( i = m c (fi)/2 is the reduced mass 
of the cc-pair, that comes upon matching between the full and effective theory decay constants 
(the coefficient Co(m re d/^) in [i^]). The exponent on the left-hand side is simply 2m c (3 c5 ■ x since 
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taking into account the binding energy of the cc-pair (see discussion below dSJ)) would amount to 
replacing 2m c — ► 2m c (l — v 2 /8), which brings corrections of ( u2 )- 

On the right-hand side of (|A6P we substitute Pjm = 2m c /3 C c = Mj/^,(3 CC . Equating both sides 
gives: 

E ( J /^fe( £ )l -f C°-px] 1°) = -*Cj/iW2M)/ W Af ^ e* (fl*). (A7) 

p x ^0 

3. Heavy-to-light current 

The Hamiltonian density 7i\y in the leading logarithm approximation is given by —C\y and its 
matrix element by (|38j) . We write it here explicitly for reference: 



(J/^ ce (e), X s \ Ctffa r b ) \B Pb ) = ^£-V cb V c * s C (to) E E E el(1 



x (x sn | g# r ) r,y„t(x) h Ph {x)\ \B Pb ) (J/^ C M\ [£ { Z P Tj vlfl p J \o P J ■ 

(A8) 



Replacing the matrix element of the cc-current with (jA7J) gives: 

71 



(J/iPf3 c M, X s \£&(», // 6 )|%) = ^|^F ( ;Co( W )[-iC^(m c /2 /1 )/^M J/f ] 



\ E {(^nl x)[ £ *(/3 cc -) • nj - ./„, A (M«); x)[e* (M ■ e +± (n)} \B Pb )} , (A9) 



x 

2 

p^O 



where 



J n ,^(r,-; X) = e mj/^+V-rn M . X[ ^ w ^) C LL {^-\ Tj hp b ]{x) (AlO) 

is the standard leading order expression for the SCET heavy-to-light current slightly modified here 
by the presence of the phase factor of the cc-state. The function C LL (n ■ V'/n) represents the 
exponential factor of ([33]) . the coefficients c[°l(r) are written explicitly (see (|3"T|) ). 

The sum over collinear labels in (|A9j) is approximated by the single term corresponding to the 
stationary point of the exponent. The label p is chosen so that the momentum components in the 
directions n and _L vanish. In other words only oscillations on the ultrasoft scale are permitted in 
the exponential factor in (|A10|) . In so doing we have to keep the residual component along fi which 
is of the same order ~ A 2 as the ultrasoft momentum of the matrix elements in (|A9[) . Expanding 
the exponent in (|A10|) in light-cone coordinates in the U-meson frame, where 2(3 b = n + n, gives: 

Tl Ti 

[Mj/4P C c -n)+n-V- m b ]- + P X + [Mj^((3 c5 ■ n) - m b ]-. (All) 
At the stationary point: 

M J/f ((3cc ■ n) + n -V - m b = and V± = 0. (A12) 
Momentum conservation in the U-meson frame reads: 

Tl Tl 

M B p b = Mj^p cS + -(E x + \p x \) + -(E x -\p x \) and Pj/ ^ = p x . (A13) 
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Using it we solve the first equation in (|A12|) for n-V: 

n-V = m b -M B + E x + \vx\ = 2Ex{l + O (A 2 )) 



(A14) 



and express the residual component along n in terms of Ej/$- Finally, taking the sum over collinear 
labels in (|A8|) amounts to replacing the effective theory current (|A10|) by the expression 

J n ,0 b (T r , X) = e *(^+|p^|-m b )K«-)£g W (0) r . y\ (A15) 

and substituting for the operator function C LL (n ■ V^/fi) its value at the stationary point 
C LL (2Ex / ' n)- Note that although V± is set to zero, the V\ is the invariant mass of the jet M x 
(see (|A3P ) and must be kept. 

Taking the square of (|A9j) and summing over polarizations of J/i/j gives: 



(J/^ cB (s), X a \C#(n, u b )\B 0b 



lG 2 F \V ch V c 



* |2 
cs I 



C (n b )C 



LL ( 2E x\ „-l 



x ') c m (,2/1 



x [M 2 B I (X sn \ J n>Pb (P R ; 0) \B Pb ) | 2 + M 2 j^\ {X sn \ J n ,p h {t-{n); 0) \B Ph ) | 2 } 



(A16) 



4. Differential decay rate 



Summing over the inclusive collinear states in the last line of (|A16|) is performed by means of 
the standard trick. The sum is equal to the imaginary part of the T-product of the effective theory 
currents (|A15P : 

- 21mT& m ) where T s {q) = -i j d 4 xe^' x ^- {B Pb \ T{J+ )A (T i; x) J„, ft (r i; 0)} \B Pb ) . 

(A17) 

Time-ordering in this expression is applied only to the fields which carry labeled momenta. The 
ultrasoft fields belong to the zero bin [42| of the sum over collinear momenta and are not affected 
by the time ordering, the corresponding operators are simply averaged over the Bp -state. 

The matrix element in (IA17I) is then reduced to a convolution of the perturbative jet-function 
and a shape function of i?-meson which is a non-perturbative object. Those functions are defined 
as: 



<0|T[W*<0) g°^ (0) ] |0> 

1 



I 



S(l 



( B l3 b \ [hp b 5 (in- D us 



d 4 k 
(2vr) 4 

l + )hp b ]\Bp b 



e - ikx J P (k + )^, 



(A18) 



Here k + = n ■ k and Z + = n ■ I. The calculation is almost identical to B — > X s j [lCj] to which the 
reader is referred (4^]. The only difference is due to Dirac structures. In our case Tj is either Pr 
or ^_(n). To average the Dirac structures over 5^-meson state we use the fact that the B-meson 
is a pseudoscalar and in HQET the spin of the heavy quark decouples from the gluon field at the 
leading order in A c /mfe-expansion. The calculation is straightforward and gives: 



and (Bg b \j 



°^(n)h°it-(n)\B Pb ) 



R i 2 ' 11 '~l J b> ^ \~ Mb I i r~\-~j i <2r—\-~j\ — p b i 2' 

Finally the differential decay rate in terms of the jet and shape functions becomes: 



(A19) 
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dE 
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47T 
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G 2 F \V cb V c 



* |2 
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C (n b )C 



LL 



(2Ex 

V V- 



\ 2a 



2 



f 2 

Jj/i/> 



M B + 2 



M B 



7T 



Im J P (k + - {E J/4 , + lpj/^1 - m b ),n) 



S(k + ,a). 



(A20) 
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The /z-dependence of the jet and shape functions cancels that one of the Wilson coefficients, so the 
differential decay rate is /i-independent. The argument of the jet function contains the invariant 
mass of the jet, V\ ~ M| given by (|A3p . The scale fi b is the matching scale of the order of the 
6-quark mass m b . 



APPENDIX B: COEFFICIENTS C^Jm b /n, r) 



i k;j ' 



Here we give the explicit expressions for the coefficients c\ l ^.-{m b /fi, r) (see (152p ). To simplify 
the output we introduce functions F\, F 2 , and F3: 



Fi(r, fi/m b ) = --Log 



m 2 . (1 - 4r 2 Y 



yLog 



TO 



vr 2 8 T . 

h -Li 2 

9 3 



1 - 4r 2 



+ -Log [l-4r 2 ] 2 , 



Fo(r) 



[1 - 4r 2 ) 



(1 - 2r 2 ) 



(Log [l - 4r 2 ] - m) , 



F 3 (r, n/m h ) 



+ 



+- 



3(1 + 4r 2 ) 
-1 + 4r 2 
17vr 2 13 



rLog 



Lis 



ml 



+ ( -Log [1 - 4r 



1 



2-4r 2 



Li 2 



2r 2 



21 6Log[2r] 13 
1 - 4r 2 3 

2r 2 



Log 



/' 



2 1 



nil 



1 - 2r 2 



3Li 2 



-1 + 4r 2 



— Lio 

3 



3i7rLog 

4r 2 
-1 + 4r 2 



(2 - 4r 2 ) 11 2 ' 
(l-4r 2 ) 2 m 2 



21 2 . o T „„ h , jn t n 0.21 12Log[2r]Log [l - 4r 2 ] 



—Log [l - 4r 2 ] + 3Log [l - 4r 2 ] Log [l - 2r 
9 3 



-1 + 4r 2 



+ 



12Log[r]Log [1 - 2r 2 ] 3 (3 + 4r 2 ) Log[2]Log [l - 2r 2 ] 6Log[r] 2 3 (3 + 4r 2 ) 



-l + 4 r 2 



-1 + 4r 2 



+ 



+ 



l-4r 2 2(-l + 4r 2 ) 



Log[2] 5 



(Bl) 



The coefficients C^, - are enumerated according to the effective theory Dirac structures j 
1,2,3. 



<$1 = \ fa*> p/m b ) - f + ^Log [1 - 4r 2 ]) , 
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80:1 
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m if, , x t7r (21 - 76r 2 + 40r 4 ) 27 - 40r 2 r , 

Cgi = ~Tr "/"") " \ (1 _, r ^ j + - 3Log[4r] 
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24r 2 (-1 + 2r 2 ) 

APPENDIX C: TRANSFORMATION OF OPERATOR BASIS 

In this section we present without derivation the relations used to calculate the Wilson coeffi- 
cients Cq LC \rrib) and C% LO (mb) given in Eq. (J3j) - The Wilson coefficients for the set of operators 

Pi = [s^P L T a c\[c^P L T a b], P 2 = [s^P L c][c^P L b}. (CI) 



have been calculated in [29j] at the NNLO in a s . We need however the Wilson coefficients in the 
different operator basis 

1 = [s ll ,P L T a b][c^P L T a c], 2 =[s^P L b][c^P L c}. (C2) 

At the LO the two bases are related by a linear transformation which follows when one applies the 
Fierz transformations in the color and spinor spaces: 

SH? !)(£)• 

At the NLO the one-loop corrections have to be taken into account and that brings compli- 
cations due to the evanescent operators. The dimensional regularization used when calculating 
Feynman integrals does not have an unambiguous definition for matrix 75, which is essentially a 
four-dimensional object. The systematic way of coping with this difficulty (see, e.g. (4p| ) consists 
in expanding the operator basis to include operators vanishing in the limit D = 4 (evanescent 
operators). In our case the ambiguity is brought by the Fiertz transformation (|C3P which is not 
defined at D 7^ 4. Working out the relation between the bases (IClj) and ()C2j) at the NLO requires 
introduction of two evanescent operators 



F 2 \0 2 V 2 i J V P 2 



Fi \ (0 1 \ ( -\ I \f P] 



3 



1 



(C4) 



and calculating one-loop diagrams analogous to those shown in Fig. [7] with the operators (|C4|) 
inserted. It suffices to extract only the UV-divergent part of the additional diagrams. 
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The Wilson coefficients of operators (|C1|) and (|C2|) are related by 



4%) r ^ 



c?V) 




(C5) 



Here C"(/i) are the Wilson coefficients for the basis (|C2p and C(/i) for the basis (IClj) . The first ma- 
trix in Eq. (|C5|) is the transposed matrix in Eq. (|C3|) corresponding to the tree-level transformation 
of the operator basis. The second matrix in (IC5h is the contribution of operators (|C4h , 

The explicit expressions for the Wilson coefficients C(fi) can be found in Eqs. (39)-(44) in (29| . 
Substituting them into Eq. (|C5j) gives the required expressions for Cq (nib) and Cg LO (rrib). 



APPENDIX D: DIRAC STRUCTURES OF THE EFFECTIVE THEORY 



In this section we present the reduction formulae necessary to project the matrix elements of 
the full theory into the matrix elements of the effective theory. The generic Dirac structure on the 
full QCD side has the form 



[sri6][cr 2 c], 

where T\ t 2 is one of the following products of Dirac matrices 

r = {i; 75; 7 M ; 7^75; 7V; 7^75; 7 M 7V; 7^7^75} . 



(Dl) 



(D2) 



Here the last four entries come from one-loop corrections to the tree-level full theory operator. On 
the effective theory side the generic matrix element is 



[fn, p rf 6 T hp J [£p cc p± r 



et c 

CC 'I Pec — pj 



(D3) 



where Tf b T and are the set of operators which form the complete set on the corresponding 
subspace of Dirac spinors. 

The Dirac structures T^| T follow by sandwiching the full theory structures of Eq. (ID2j) between 
the projectors Pp cE and P-p cS 



(D4) 



Working out the Dirac structures Tf b T requires more effort because now the original Dirac structure 
from the list (1D2|) is sandwiched between the projectors on different subspaces: Pn and Pp b . The 
most straightforward way to work out the necessary relations is using Eq. (17) from [4lJ that reads 



r 



ET 



2 



4 



^PnVPp 



4 m 

|75Tr 



^PnTPp b 



+ 7x Tr 



luPnTPp 



(D5) 



n' 



n 1 



'1 

2 • 



Here 7^ = 7^ 

The projection formulae are considerably simplified in the restframe of the 6-quark (6RF) where 
2(5b = n + n. The four-momentum of the s-quark in the 6RF is specified by its energy: = EgU^. 
To simplify the results further we use the relations 



7± 



• I1U 



-2e^_, 



7j_ + l£ ± IvL 



-2et 



(D6) 



where e± are the polarization vectors of s-quark with ± corresponding to the spin direction 



along/opposite to the quark momentum. Here £j_ 
ified explicitly in the 6RF where = (1,0,0,1) 



£^ ur iP P b n p . These relations can be ver- 



= 2 1//2 (0, 1, ±i, 0). Equations below summarize the results. 



(1,0,0,-1) 



(1,0,0,0), and 
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On the cc-side the following relations hold: 

® 7 &y]p fl ®\{%<nl - fei) + ® l -e^P c5p7v± . (D7) 
Here 7^ = 7^ — and 7^7^ = ^[7^,7"]. On the s6-side we use relations 



7 [^ 7 Hp i(g) ( ra ^ _ n^ 6 M )Pc ® +ief Pi ® +(n^ - nV 
7 [/ V ] ^® - W " ""flftPjl ® -ieTPfl ® +K4 - nVM^_ <g> . (D8) 



Using Eqs. (|D7[) and (|D8j) it is straightforward to verify (|18p . Finally, there is one more relation 
holding at the leading order in f-expansion that allowed us to reduce the number of the effective 
theory Dirac structures when processing the output of the full QCD calculation: 



E, 



®«f±. (D9) 

It follows from conservation of momentum, m&/3& = 2m c /3 cE + E s n, and observation that on the 
cc-side only the components perpendicular to (5 C c contribute, which is obvious from (|D7D . 
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